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Abstract 

It is shown that the dynamics of cosmologies sourced by a mixture of per¬ 
fect fluids and self-interacting scalar fields are described by the non-linear, 
Ermakov-Pinney equation. The general solution of this equation can be ex¬ 
pressed in terms of particular solutions to a related, linear differential equation. 
This characteristic is employed to derive exact cosmologies in the inflationary 
and quintessential scenarios. The relevance of the Ermakov-Pinney equation 
to the braneworld scenario is discussed. 
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1. Introduction: Observations of the cosmic microwave background (CMB) power 
spectrum [|I| now provide strong support for the inflationary scenario p| (for recent 
reviews, see, e.g., |0). The simplest mechanism for inflation utilises the potential 
energy associated with the self-interaction of a scalar inflaton held to drive the accel¬ 
erated expansion. High redshift observations of type la supernovae suggest that the 
universe is experiencing another phase of accelerated expansion at the present epoch 
This, combined with the CMB data, presents a picture of the universe dominated 
by a dark energy component ii- One possible source of this dark energy is a scalar 
quintessence held that interacts with baryonic and non-baryonic matter in such a 
way that its potential energy is currently dominating the cosmic dynamics 0. The 
observations favour a model of structure formation where 70% of the energy density 
of the universe is presently in the form of quintessence |p. The remaining fraction of 
the energy density is comprised of visible and cold dark matter which collectively act 
as a pressureless perfect huid. 

The ekpyrotic scenario has recently been proposed as an alternative to the stan¬ 
dard inhationary cosmology [^. In this scenario the big bang is interpreted as the 
collision of two domain walls or branes travelling through a hfth dimension. Be¬ 
fore the collision, the ehective dynamics on the four dimensional branes is described 
by Einstein’s gravity minimally coupled to a self-interacting scalar held. The held 
parametrizes the separation between the branes and at early times slowly rolls down 
a negative potential. This results in an accelerated collapse of the universe and since 
the held is minimally coupled, its energy density is related to the Hubble parameter 
by the standard Friedmann equation. Thus, self-interacting scalar helds play a cen¬ 
tral role in modern cosmology and in view of the above developments, it is important 
to investigate cosmologies that contain both a scalar held and a perfect huid. 

In this paper, we develop an analytical approach to models of this type by express¬ 
ing the cosmological held equations in terms of an Ermakov system [^, |^, |^. In 
general, an Ermakov system is a pair of coupled, second-order, non-linear ordinary 
diherential equations (ODEs) [|TI| and such systems often arise in studies of nonlinear 
1^, nonlinear elasticity ||^, molecular structures 


PI 


quantum held theory 
(For further references, see, e.g.. 


optics 

in curved spaces |]^ and quantum cosmology |l6 
Refs. [0). 

In the one-dimensional case, the two equations decouple and the system reduces 
to a single equation known as the Ermakov-Pinney equation |^, |^. This is given 

b>f| 

^ + ^=^ 3 . ( 1 ) 

where Q is an arbitrary function of r and A is a constant. Although Eq. (J^) is a 
non-linear ODE, it exhibits a remarkable superposition property [IT, that implies 
that its general solution can be expressed directly in terms of particular solutions to 
the related linear, second-order ODE, where A = 0 
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'Eq. § is sometimes referred to as the Milne-Pinney equation 


1 










































2. The Ermakov-Pinney Equation in Cosmology. We begin by deriving the 
Ermakov-Pinney equation in a cosmological context. The held equations for a spa¬ 
tially hat, Friedmann-Robertson-Walker (FRW) universe with a scalar held and per¬ 
fect huid matter source are given by 

= + + (2) 

0 + 3ih0+^ = O, (3) 


where = 0^/2 -|- R(0) is the energy density of the scalar held with potential R(0), 
Pmat = Da~^ is the energy density of the barotropic perfect huid with equation of 
state pmat = [(^ —3)/3]pmat, & dot denotes diherentiation with respect to cosmic time, 
t, H = a/a represents the Hubble parameter, a is the scale factor of the universe, D 
is an arbitrary, positive constant, 0 < n < 6, hC = and mp is the Planck mass. 

Eq. represents the conservation of the energy-momentum of the scalar held and 
can be expressed in the form = —?)H(jP. 

One method of reducing this second-order, non-linear pair of ODEs (^-(^ is 
to write them as a hrst-order system. This can be done if the scalar held varies 
monotonically with cosmic time. It then follows that Eq. can be expressed as 

^ = -'iHi' ( 4 ) 

and the Friedmann equation then takes the form 

where y = a"". The perfect huid source leads to a non-trivial right-hand side in Eq. 
(|). The general solution to Eq. (^) can be expressed in terms of quadratures if the 
functional form of p<^(0) is known: 
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a(0) = exp 




X <( n — ntCD [ dcf) 
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Ijn 


( 6 ) 


where H is a constant of integration. The ehects of the huid source are contained in 
the second term. 

The disadvantage of this change of variables is that it is not valid if the scalar held 
exhibits oscillatory behaviour at any time. We therefore consider a diherent approach 
by diherentiating Eq. and substituting Eq. (|^). This implies that 
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a a? 



(7) 
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We now define an effective scale factor, b: 


and a new time parameter, r: 


= 62/- 


d _ 

dt dr 


( 8 ) 

( 9 ) 


and since b is positive-definite, r is a monotonic increasing function of cosmic time, t. 
Eq. ( 0 ) then transforms into the second-order, non-linear Ermakov-Pinney equation 


d% ^ riK^ (dcj)'' 
dr , 


b = 


4 

and comparison with Eq. ([^) then implies that 

.2 / . - \ 2 
V dr 


Dn^K^ 1 


( 10 ) 


Q = 


UK" I d(j) 


A = 


Dn‘^K‘^ 

1^' 


( 11 ) 


The change of variables (|^) and (0) is interesting because Pinney |T^ has shown 
that two linearly independent solutions, ( 6 i(r), 62 ('r)), of the time-dependent oscilla¬ 
tor equation 

— -h Q(r )6 = 0 ( 12 ) 

can be combined to give the general solution to the Ermakov-Pinney equation (0): 

1/2 


bp = 


Abi + Bb ^ + 2^6162 


(13) 


where {A,B,C} are constants satisfying the constraint equation 

AB-C^ = 

and the Wronskian 

W = bi — 
dr 

is a constant due to Abel’s theorem. It can be verified by direct differentiation that 
Eq. (^) does indeed satisfy Eq. (i- An important property of Eq. ([T^) is that if 
a particular solution, 61 (r), is known, a second solution can be written in terms of a 
quadrature: 

r'’' Ht 

h(T) = Wbi(T) 


A 

(14) 

dbi 

(15) 




r 


(16) 


Thus, the general solution to the Ermakov-Pinney equation (pT]) is found in terms of 
a particular solution to Eq. (E3)- 
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The perfect fluid component {D ^ 0) results in the non-linear sector on the right- 


hand side of Eq. ([10|) and the effects of the scalar held are parametrized entirely in 
terms of the function, Q, dehned in Eq. O- We refer to this function as the 
‘oscillator potential’. It is worth noting that Eq. (^) is independent of the specihc 
functional form of the scalar held interaction potential, E( 0 ). 

For the pure scalar held model {D = 0), Eq. (|T^ may also be interpreted as a 
classical analogue of the zero-energy Schrodinger equation, where the scale factor, 
b, plays the role of the ‘wavefunction’ and the ‘wave vector’, Q(r), is determined by 
the kinetic energy of the scalar held. An equation of this form has arisen previously 
in studies of FRW cosmologies containing only a perfect huid source with a constant 
equation of state oj = Pmat/Pmat lUl and for the special case of a self-interacting 
scalar held that satishes this condition [^. This analogy may be extended further 


by considering an alternative reparametrization, where a new scale factor and time 
variable are dehned by c = and cr = f dtc, respectively. It follows that Eq. (0) 

may then be expressed as a one-dimensional Schrodinger equation: 


cPc 


+ [E - P{a)] c = 0, 


(17) 


where E = —{'n?K^D)/12 represents the total energy of the system, and is determined 
by the momentum of the perfect huid, and the potential energy is given by P = 
{uk'^/ 4){d(j)/da)‘^ in terms of the scalar held. 

There is a close relationship between the Schrodinger equation ( 0 ) and the 
Ermakov-Pinney equation (|^), as discussed by Milne |]^, who developed a method 
for solving the former equation in full generality in terms of a particular solution of 
the latter, in such a way that the oscillatory behaviour of the wavefunction is mani¬ 
fest. In this approach, the general solution to Eq. (JT^), or equivalently Eq. (0) , is 
expressed in the form 

bn = cos (yX6{T) + , (18) 

where {U,e} are arbitrary constants and {'0 (t), 6 '(r)} are functions that need to be 
determined. Substituting this ansatz into the homogeneous equation (0) reveals that 
Eq. ([ISD represents the general solution if these functions satisfy the coupled system 
of ode’s: 


+ ^ 


dr^ 


de 

dr 


1 


(19) 

( 20 ) 


and it follows, therefore, that the general solution to the Schrodinger equation (|^ 
can be deduced, at least with respect to the quadrature, 0 = f'^ df'ijj~‘^ {f) , if a solution 
to the Ermakov-Pinney equation O is known. 

Furthermore, two linearly independent parametric solutions to Eq. © are given 
by bi = ^jJ cos (y/XO) and 62 = '^sin (a/A^), where the Wronskian is = A, 9{t) 
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satisfies Eq. (^OD and t/’(r) is a particular solution to Eq. (^). Substituting these 
two solutions into Eq. (|T^ gives 


bp — 


ip [hlcos^ (VA0) + Esin^ (yx^) + {AB - (2^6) 


1/2 


( 21 ) 


where we have employed Eq. (W). This implies that the general solution to the 
Ermakov-Pinney equation (|^ can be determined in terms of a particular solution, 
pj^r), to this equation 
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To summarize thus far, the dynamics of a pure scalar held cosmology is deter¬ 
mined by a one-dimensional oscillator equation with a time-dependent frequency. 
Depending on the parametrization chosen, the inclusion of a perfect huid source re¬ 
sults in an effective constant shift in the total energy of the system, as in Eq. ([T7|) , 
or equivalently introduces a non-linear term as is the case in Eq. m- In particular, 
the problem of solving Eq. (|T0|) for a coupled perfect huid and scalar held cosmology 
has been reduced to solving Eq. ([l^ for a given functional form of Q{t). This is 


interesting because this latter equation has been studied extensively in the literature. 
For example, algorithms that extend the Wentzel, Kramers and Brilluen (WKB) tech¬ 
nique have recently been developed to hnd power series solutions to Eq. (^) [p^ . 
For related work, see also Ref. |2^ . 

3. Algorithms for Solving the Field Equations: We now proceed to outline an 


algorithm for deriving exact cosmological models when both a scalar held and perfect 
huid are dynamically signihcant. Attention to date has primarily focused on pure 
scalar held models 
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28|. Some specihc models with a perfect huid were 


also recently analyzed |^. /^From a particle physics perspective the potential of the 
scalar held is the fundamental quantity that dehnes the model. In general, however, 
Eqs. (H) and are very difficult, if not impossible, to solve unless the potential has 
a specihc form, as in the exponential case [^. There are several possible routes to 
take when solving the Ermakov-Pinney equation (|I[). Given its properties as discussed 
above, one possibility is to specify the time dependence of the scalar held, (p = 
since the functional form of Q = Q{t) is unaltered when deriving the solution (|^ 
from the corresponding solutions to the homogeneous equation (|^). 

/,From the cosmological point of view, however, it is the time-dependence of the 
scale factor that is important. In particular, the ehective equation of state and the 
potential of the scalar held can in principle be determined directly from high redshift 
observations if the dependence of the scale factor on redshift is known to a sufficiently 
high accuracy [^, |^. This question has recently attracted renewed interest in light 
of the proposed Supernova/Acceleration Probe (SNAP) ^ . 


Instead of choosing the potential, therefore, an alternative approach is to invert 
the problem by hrst specifying the dependence of the scale factor on cosmic time 
mM- We may begin with an ansatz for the scale factor, ai{t), or equivalently, the 
rescaled function, bi(t), and use Eq. (^) to hnd Q{t). Integrating this function then 
yields the time-dependence of the scalar held. The general solution to the Ermakov- 
Pinney equation can be deduced immediately from Eqs. (pK]) and ([T6|). The scale 
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factor for the Ermakov-Pinney equation in parametric form then follows from Eq. 
(I). Its dependence on cosmic time may be evaluated by integrating Eq. (H) and 
inverting the result. Finally, the scalar held potential is reconstructed directly from 
the Friedmann equation (H), rewritten in terms of the new variables (H) and 


V{t) 


12 _ 1,2 f^y _ £_ 


( 22 ) 


The form of E((/)) then follows by inverting 0(r) and substituting r^cf)) into Eq. (p^). 
The various possible algorithms are illustrated in Figure 1. 



Figure 1: Illustrating the possible routes for finding the scalar field potential in a pure 
scalar field system and in a system comprised of both a scalar field and a perfect fluid. The 
cosmological Ermakov-Pinney equation ([To|) plays a central role. The two potentials can 
be found from a given functional form of the oscillator potential, Q(r), or from the scale 
factors, a{t). In principle, it is consistent to begin at any point on this Figure. From a 
cosmological perspective, the scale factor is the primary function of interest, and specifying 
this parameter from the outset is physically well motivated. 

4- Exact Cosmologies Sourced by Scalar Fields and Perfect Fluids: We now illus¬ 
trate this procedure with a number of examples. The simplest case is that of the scale 
factor as an exponential function of time ai{t) = where the power has been cho¬ 
sen so that &i(r) = r and therefore fe 2 ('r) = 1- Such a choice produces a static scalar 
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field, dcp/dr = 0, corresponding to a pure cosmological constant. The solution for 
a cosmological constant with a perfect fluid is already known [^, but this example 
illustrates the generality of the technique. Thus, from Eqs. (|]) and (|^, the general 
solution for the scale factor when a perfect fluid is present is given by 


ap = 


{At^ 


+ 2Ct + B 


l/n 


(23) 


where AB — C^ = — Dri^n^/12. Since the scalar held remains independent of time, the 
potential is also a constant, V = 12A/(K^n^), from Eq. (|^) . Finally, the dependence 
of the scale factor on cosmic time is deduced by using Eq. and integrating Eq. 

(I): 

l/n 

sinh^/" (^^/At) . (24) 

A more general ansatz for the scale factor is given by 



a±(t) = Gf''±, (25) 

where q± are constants and G = [2/(2 + ng±)]^±. Using the relationship between a 
and b given in Eq. (^, it is possible to express the two solutions to the homogeneous 
equation as 

= (26) 

where p± = nq±/{2 + nq±). Solving the homogeneous equation Eq. (^) then gives 
the form of the scalar held 4>{t), 


0 = Flnr, (27) 

where the constant F is related to the power p± via the quadratic equation 

P± = --• (28) 


The Wronskian is = l—nn^F^ and, for consistency, we assume that nn^F^ < 1 
in what follows. The general solution to Eq. (|T0|) is then found by substituting Eq. 
(p6|) into Eq. (O) and is 


CLp — (Ar2p+ + Bt^p- + 2Ct) 


l/n 


(29) 


The analytical form of the potential deduced from Eq. (^^ is a combination of 
exponential functions of the scalar held. The simplest case arises for A = B = 0. 
Eqs. O and ([T^) then imply a relationship between the constants: 


D = 


12C^ 


(^1 — nn'^F'^^ 


(30) 
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and consequently, the scale factor, the scalar held and its potential are given by 


ap{t) = 

(p{t) = F\n +2F\nt 


V{4>) = 


CF^{6 - n) 


n 



(31) 


The solution (|^) describes a cosmology where the energy density of the scalar 
held tracks that of the perfect huid in a way that leaves the dependence of the scale 
factor on cosmic time unaltered from the pure perfect huid model 0. The self¬ 
interaction coupling, F, of the scalar held is constrained only by the condition that 
the huid integration constant, D, be real. 

Eq. (^) represents a past or future attractor for the more general classes of 


solutions given by Eq. (|29|) . For example, if A ^ 0 and B = 0 the potential (shown 
in Figure 2 for the case of a pressureless huid) is given by 


-7 


12 




,-20/F 


D 


(32) 


where 


Y = + 2Ce^/^ 

Z = + C. (33) 


In the early-time limit (r —0), Eq. 
hand, the potential reduces to 


(|3^ asymptotes to Eq. (^). On the other 


Eoo 


A 


24pl 





exp 



(34) 


at late times (r —>• cxo), since > 1/2. In this limit, the scale factor asymptotes to 
a cx where p = 2p^/{np_) and it is straightforward to verify that this power law 
behaviour is the late-time attractor for a universe dominated by a single scalar held 
with an exponential potential, V oc exp {^\/2K(f )/. Since > p_, the power of the 
expansion in this limit, p, exceeds that of the pure perfect huid model (a cx t^/”) and 
it is greater than unity if p+ > nnF/ -\/8. In this region of parameter space, therefore, 
the solution given by Eqs. (p^ ) and ( |3^ describes a universe that behaves as a perfect 
huid cosmology at early times. However, the perfect huid becomes negligible as the 
expansion proceeds and the universe subsequently enters a phase of power law, inha¬ 
tionary expansion driven by the exponential potential of a scalar held. For n = 3, the 
huid is pressureless and current observations indicate that the universe is undergoing 
such a transition at the present epoch [^ |^, |[. Thus, it has been demonstrated in 
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Figure 2: The potential (^) for ^ = 1, S = 0, C = Vs, D = 1, F = 1/2, k = 1 and n = 3. 
The potential asymptotes to an exponential form at high and low energies. 


the example outlined here that the algorithm for solving the cosmological held equa¬ 
tions (H) and (H) using the Ermakov-Pinney equation (|l^) has generated a new exact 
solution describing the quintessence scenario. Moreover, exponential potentials often 
arise in superstring-inspired models through non-perturbative effects . 


Another class of solvable models arises when the scale factor is a trigonometric 
function of r. The two solutions to Eq. (^) are bi = sin yr and 62 = cos yr, where y 
is a constant and the Wronskian is PF = —y. Given this ansatz for the scale factor, 
the scalar held varies linearly with r, such that 
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^/nK 


T. 


The solution to Eq. (pT]) is then given by 


Cbp — 


A sin^ 'jT + B cos^ yr -|- G sin 2yr 


l/n 


(36) 


(36) 


In general, the potential is a complicated expression of trigonometric functions. How¬ 
ever, for the particular choice A = B = 0, Eq. (^) takes the simple form 


V = 


2Cj^ 



(37) 
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and it follows from Eq. (p6D that V oc ap. Consequently, the potential energy of the 
scalar held is inversely proportional to the energy density of the perfect huid. 

So far in this section it has been demonstrated that once a solution to the pure 
scalar held system given in Eq. (|T^ is known, a solution to the Ermakov-Pinney 


equation (0 describing a scalar held and perfect huid cosmology can be found. 
Before concluding this section, we employ the correspondence summarized in Eqs. 
(P!^ ) -(PI1| ) to show the reverse procedure, i.e., to derive a pure scalar held cosmology 
from a known solution to the Ermakov-Pinney equation (^). We invoke the ansatz 
= Ae^* which in terms of the rescaled time variable t is 'ip = At, where A is a 
constant. Substituting this ansatz into Eq. (|^ and solving Eq. (^) for the scalar 
held gives 


= T 


4A 


TIK^ 


1/2 


AV’ 


(38) 


where the integration constant is specihed to be zero without loss of generality. In¬ 
tegrating Eq. (1^) implies that 6 = — l/(A^r) and after specifying U = 1 and e = 0 
for simplicity, the scale factor of the universe and the potential for such a pure scalar 
held model are given by 


a = 


At cos 


A^r 


P = 


2A2 r 




6 (cos(p -t- (psin(p)" 


rup^ cos^ (p 


(39) 


(40) 


where = A/nK0/2. In the small held limit, (p -C 1, the potential (^) approximates 
to an inverted harmonic oscillator and for sufficiently small values of the scalar held, 

m 


the universe undergoes inhation 


Potentials of this form generically produce a 


spectrum of density perturbations in the post-inhationary universe that increases 
in amplitude on larger scales without producing a detectable signal of gravitational 
Observations of CMB anisotropies presently favour such a spectrum 


waves 
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5. Braneworld Scenarios: It is also of interest to consider whether the approach 
outlined above can be applied to other cosmological scenarios. Considerable attention 
has been focused recently on the possibility that the observable universe may be 
viewed as a domain wall or brane that is moving along a timelike geodesic in a 
hve-dimensional, static, bulk spacetime |40, 41, 42, where the brane equations of 


motion are determined by the Israel junction conditions [44]. An observer conhned to 


the domain wall interprets such motion in terms of cosmic expansion or contraction. 
The junction conditions relate the second fundamental form (extrinsic curvature) of 
the induced metric on the brane to the energy-momentum tensor of matter conhned 
on the brane. For the case of a spatially hat brane moving through a Schwarzschild- 
anti de Sitter space, the ehective Friedmann equation on the brane is, after some 

I m m m 


simplifying assumptions 




= -p [I + 


given by 
P 


2A 


+ 


pn 


(41) 


10 














































where A is the brane tension, is related to the mass of the black hole in the bulk 
and n = 4. The standard, linear dependence of the Friedmann equation on the energy 
density, p, is recovered in the low-energy regime, p -C A. Further modihcations to 


the Friedmann equation (^I]) may also arise in more general braneworld settings. For 
example, a term of the form Eg = where pg < 0, arises if the black hole carries 

an electric charge [^, . 


For this model, the energy-momentum tensor of matter on the brane is covariantly 
conserved, and consequently, Eq. (|^) must also be satished for the case of a single, 
self-interacting scalar held. Eqs. (|^) and (|4T|) are therefore relevant to a number 


of braneworld scenarios. By differentiating Eq. (|4TD with respect to cosmic time. 


substituting for the scalar held from Eq. (y), and using the change of variables (§) 
and (^, it can be shown that the dynamics of this system can be written in the 
Ermakov-Finney form 

d% ~ _ 

—— -\- Qb — 




where 




UK 


1 + 


n fin 

(42) 


(43) 


In Eq. the ehects that lead to the quadratic dependence of the Friedmann 

equation ( ^ID on the energy density of the scalar held modify only the oscillator 
potential (|4^) , whereas the non-linear term on the right-hand side arises through the 
bulk black hole contribution. Thus, in principle, a similar approach to that outlined 
above in Section 3 may be developed for analytically investigating the importance of 
these terms on braneworld cosmologies. For example, if a given solution is known 
for the case where there is no black hole in the bulk |^, it may serve as a seed 
for generating a corresponding cosmology when the black hole is present. This is 
interesting since it has recently been shown that such terms can signihcantly alter 
the qualitative dynamics of braneworld cosmologies . 


6. Summary. In this paper it has been demonstrated that after a suitable re- 
dehnition of variables, Eqs. (@)-(|3) can be related to the simplest Ermakov system, 
thus reducing the problem of solving these equations to solving a single second-order, 
linear ODE. Moreover, the nature of the Ermakov-Pinney equation implies that there 
exists a correspondence between a spatially flat FRW universe containing a scalar held 
and a cosmology containing both a scalar held and a perfect huid. The functional 
form of the scalar held, 0(r), is identical for the two models, but the potentials, V{(f)), 
are diherent in each case. We have proposed an algorithm for hnding exact solutions 
for a scalar held and perfect huid model by employing a pure scalar held cosmology 
as a seed. The algorithm is founded on the property that a particular solution to the 
homogeneous equation (|T^ leads to the general solution of Eq. (Q). It has also been 
shown that a model containing both a scalar held and perfect huid can act as a seed 
for hnding more general pure scalar held cosmologies, since a particular solution to 
Eq. (2) allows the general solution of Eq. (|^ to be deduced |]^. Exact solutions 
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are important because they allow one to gain insight into the generic nature of cos¬ 
mologies of this type and they also provide a framework for classifying the different 
behaviours that may arise. 

The emphasis of the present work has focused on the effect of introducing a perfect 
fluid source. However, the approach we have developed may also be employed to 
consider pure scalar held cosmologies in more general spacetimes. In particular, the 
last term on the right-hand side of Eq. may be interpreted as arising from the 
spatial curvature of the universe if n = 2, where D > 0 (D < 0) for negative (positive) 
curvature. On the other hand, the shear of the spatially hat, anisotropic Bianchi type 
I metric leads to a term where n = 6. 

Finally, it would be of interest to investigate whether other classes of scalar held 
cosmologies can be analyzed in terms of an Ermakov system. Scalar-tensor cosmolo¬ 
gies provide one such example, where the scalar held plays the role of the gravitational 
coupling IQ . Many higher-order and higher-dimensional theories of gravity, includ¬ 
ing the string ehective actions, may be expressed in a scalar-tensor form |^. Since 
these theories are conformally equivalent to Einstein gravity that is minimally cou¬ 
pled to a scalar held, it is to be expected that a similar approach can be employed. 
A related, but distinct, class of theories has recently been developed to allow for a 
possible variation in the hne-structure constant |Q] . Observational evidence for such 
an ehect has been growing Q] and this behaviour can again be parametrized in terms 
of a scalar held that is coupled to the matter helds in an appropriate fashion. 

In conclusion, therefore, Ermakov systems have applications in many branches of 
mathematics and physics. In this paper, we have found that the simplest such system 
is relevant to scalar held cosmological models that are favoured by astrophysical 
observations. 
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